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Abstract 



We study the nature of the multicritical point in the three-dimensional 
0(3)©0(2) symmetric Landau-Ginzburg- Wilson theory, which describes the 
competition of two order parameters that are 0(3) and 0(2) symmetric, 
respectively. This study is relevant for the SO(5) theory of high-T c su- 
perconductors, which predicts the existence of a multicritical point in the 
temperature-doping phase diagram, where the antiferromagnetic and super- 
conducting transition lines meet. 

We investigate whether the 0(3)00(2) symmetry gets effectively enlarged 
to 0(5) approaching the multicritical point. For this purpose, we study the 
stability of the 0(5) fixed point. By means of a Monte Carlo simulation, we 
show that the 0(5) fixed point is unstable with respect to the spin-4 quartic 
perturbation with the crossover exponent ^4 = 0.180(15), in substantial 
agreement with recent field-theoretical results. This estimate is much larger 
than the one-loop e-expansion estimate ^4 = 1/26, which has often been 
used in the literature to discuss the multicritical behavior within the SO (5) 
theory. Therefore, no symmetry enlargement is generically expected at the 
multicritical transition. 

We also perform a five-loop field-theoretical analysis 
of the renormalization-group flow. It shows that bicritical systems are not 
in the attraction domain of the stable decoupled fixed point. Thus, in these 
systems — high-T c cuprates should belong to this class — the multicritical point 
corresponds to a first-order transition. 

PACS Numbers: 74.20.-z, 74.25.Dw, 64.60.Kw, 05.70.Jk 
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I. INTRODUCTION 



Multicritical phenomena arise from the competition of distinct types of order. 1-4 More 
specifically, a multicritical point (MCP) is observed at the intersection of two critical lines 
characterized by different order parameters. MCPs are expected in the temperature-doping 
phase diagram of high-T c superconductors, since, at low temperatures, these materials ex- 
hibit both superconductivity and antiferromagnetism depending on doping. The SO (5) 
theory 5 ' 6 attempts to provide a unified description of these two phenomena by introducing a 
three-component antiferromagnetic order parameter and a ci-wave superconducting complex 
order parameter. This theory predicts a MCP in the temperature-doping phase diagram 
where the two order parameters become both critical. 

Neglecting the fluctuations of the magnetic field and the quenched randomness intro- 
duced by doping — see, e.g., Ref. 7 for a critical discussion of this point — the critical behavior 
at the MCP is determined by the 0(3)©0(2)-symmetric Landau-Ginzburg- Wilson (LGW) 
Hamiltonian, 8-14 given by 

U = J d d x[\{d,ct>) 2 + \{d^f + \r^ 2 + \r^ 2 + > o (0 2 ) 2 + yM^ 2 ) 2 + > O 0V], (1-1) 

where and tp are respectively three- and two-component real fields, associated with the 
antiferromagnetic and superconducting order parameters, respectively. The critical behavior 
at the MCP is determined by the stable fixed point (FP) of the renormalization-group (RG) 
flow when both and are tuned to their critical value. An interesting possibility, put 
forward in Refs. 5, 6, 8, 9, is that the 0(3)©0(2) symmetry gets effectively enlarged to 0(5) 
when approaching the MCP. This requires the stability of the 0(5) FP in the theory (1.1). 

Cuprates have a pronounced two-dimensional layer structure with relatively weak cou- 
plings between adjacent Cu02 planes, so that two- and three-dimensional models should 
represent two extreme conditions for the possible range of properties of real high-T c super- 
conductors, see, e.g., Ref. 14. In this paper we investigate the stability properties of the 
0(5) FP in the three-dimensional (3-d) multicritical theory (1.1). 

The phase diagram of the model with Hamiltonian (1.1) has been investigated within the 
mean- field approximation in Ref. 1 (see also Ref. 13). This analysis predicts the existence 
of a bicritical or tetracritical point. The nature of the MCP depends on the sign of the 
quantity A = u v — 9-w 2 ,: If A > the MCP is tetracritical as in Fig. 1, while for A < 
it is bicritical, as in Fig. 2. 

The nature of the MCP in the 0(3) ©0(2) LGW theory has been extensively studied. 
In Refs. 5, 8, 9, it was speculated that the MCP is 0(5)-symmetric. In three dimensions, 
this picture is supported by Monte Carlo (MC) simulations of a five- component 0(3) ©0(2)- 
symmetric spin model 9 and by a MC study of the quantum projected S0(5) model. 15 These 
numerical studies showed that, within the parameter ranges considered, the scaling behavior 
at the MCP is consistent with an 0(5)-symmetric critical behavior. On the other hand, the 
one-loop e-expansion results of Ref. 3 show that the 0(5) FP is unstable. Within this 
approximation, the stable FP is the biconal FP, which has only 0(3) ©0(2) symmetry. 

The 3-d RG flow of 0(ni)©0(n 2 ) symmetric LGW theories for generic n\ and n 2 was in- 
vestigated within the e expansion to five loops. 4 This analysis showed that the 0(5) FP and 
the biconal FP are both unstable. The stable FP is the decoupled FP (DFP) corresponding 
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FIG. 1. Phase diagram with a tetracritical point. Here, T is the temperature and g a second 
relevant parameter. 

to a multicritical behavior in which the two order parameters and ijj are effectively uncou- 
pled. The stability of the DFP can be proved by nonperturbative arguments. 16,7 Indeed, the 
RG dimension y w of the operator <f> 2 i/j 2 that couples the two order parameters is given by 

y w = — + — - 3, (1.2) 

^0(3) VO(2) 

where fo(3) an d ^0(2) are the critical exponents of the 0(3) and 0(2) universality classes. 
Using 17 ' 18 i/ 0( 3) = 0.7112(5) and i/ 0(2) = 0.67155(27), one finds y w = -0.1048(12). Note also 
that the MCP described by the DFP is always tetracritical, as in Fig. 1. No other stable FP 
is found within the e expansion for the 0(3) ©0(2) theory (1.1). 4,19 According to these RG 
results, the asymptotic approach to the MCP within the SO (5) theory must be characterized 
by a decoupled tetracritical behavior or by a first-order transition if the system is outside 
the attraction domain of the DFP. 

The analysis of Ref. 4 indicates that the 0(5) symmetry can be asymptotically realized 
only by tuning an additional relevant parameter, beside the double tuning required to ap- 
proach the MCP. Of course, if the system parameters are close to those needed to obtain an 
0(5) MCP, one observes a crossover from the 0(5) behavior to the asymptotic one, which 
could be a decoupled critical behavior or a first-order transition. This crossover can be 
described in terms of universal scaling functions determined by the 0(5) FP. For instance, if 
To (5) is the critical temperature of the 0(5) MCP, g 2 and g 4 are scaling fields (functions of 
the system parameters) associated with the two relevant perturbations of the 0(5) FP such 
that g 2 = #4 = at the 0(5) MCP, then the singular part of the free energy can be written 
as 

F sins = t d »f(g 2 t-^,g 4 r^ 4 ), (1.3) 

where t — (T — T ( 5 ))/T (5) is the reduced temperature. Here v is the 0(5) correlation- 
length exponent, 02,2 > and 4i4 > are the crossover exponents associated with the two 
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FIG. 2. Phase diagram with a bicritical point. Here, T is the temperature and g a second 
relevant parameter. The thick line ("flop line") corresponds to a first-order transition. 

relevant perturbations #2 and #4, which are determined by the 0(5) FP itself. 

It has been noted 16 ' 15 ' 14 that the crossover exponent 4i4 associated with this additional 
instability of the 0(5) FP is rather small. This may partially explain the apparent 0(5) 
multicritical behavior observed in Refs. 9, 15: The MC simulations are still probing a region 
in which the crossover towards the asymptotic critical behavior is so slow to be undetectable. 
Within the e expansion, one finds 4i4 = + 0(e 2 ). Using the value suggested by setting 
e = 1, i.e. 4j4 ~ 1/26 = 0.038..., the authors of Refs. 15, 14 argued that any significant 
deviation away from the 0(5)-symmetric FP can be observed in experiments and realistic MC 
simulations only when the reduced temperature is t cross ~ 10" 10 , making the departure away 
from the 0(5) symmetric point practically unobservable. On the other hand, the analysis 
of the perturbative expansions in two different schemes, the 3-d massive zero-momentum 
scheme (six loops) and the e expansion (five loops), gives a much larger estimate of 4j4 , 4 
4j4 Ki 0.15, which makes the statement of Refs. 15, 14 on the effective relevance of the 0(5) 
FP very questionable. 4 

In this paper we return to the issue of the stability of the 0(5) FP. We compute the 
crossover exponent 4>4 using lattice techniques based on MC simulations. We obtain 4j4 = 
0.180(15). This result is slightly larger than the field-theoretical (FT) ones reported in Ref. 4 
and fully confirms the fact that the naive application of the one-loop e-expansion result gives 
an estimate that is unrealistically small. Moreover, we perform a detailed FT analysis of 
the RG flow in the 3-d scheme known as minimal-subtraction scheme without e expansion. 20 
We find that bicritical systems that have Ao < never flow towards the DFP. Therefore, 
these systems, which are those of interest for high-T c superconducting materials according 
to the analysis of Ref. 11, should have a first-order MCP. 

The paper is organized as follows. In Sec. II we discuss the stability of the 0(5) FP and 
present a MC determination of the crossover exponent 4)4 . In Sec. Ill we discuss the RG 
flow and determine the attraction domain of the DFP. In Sec. IV we report our conclusions. 
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The five-loop perturbative series used in Sec. Ill to determine the RG flow are reported in 
the appendix. 



II. STABILITY OF THE 0(5) FIXED POINT 



A. General considerations 



In order to determine the stability properties of the 0(5) FP, we must determine the 
RG dimensions of the perturbations present in the Hamiltonian (1.1) that break the 0(5) 
symmetry down to 0(3) ©0(2). They can be expressed in terms of homogeneous polynomials 
of the fields P m i, where m is the power of the fields and I the spin of the representation of 
the 0(5) group. The classification in terms of spin values is particularly convenient, since 
polynomials with different spin do not mix under RG transformations and the RG dimension 
y m i does not depend on the particular component of the spin-/ representation. We refer to 
Ref. 4 for details. For m = 2 (resp. 4), the only possible values of I are Z = 0,2 (resp. 
I = 0,2,4). Explicitly, we define 

P 2 ,O = |(0 2 +^ 2 ), 

^ 2 = W - 

Pl,0 — °2,0> 

P 4 2 = P2 0P2 2 1 

P M = §0V-^(^) 2 -i(0 2 ) 2 . (2.1) 

Note that in general all these operators renormalize multiplicatively except P 4i 2 that may 
mix with the lower-dimensional operator P 2j2 . 

Hamiltonian (1.1) can be rewritten in the form 

H = j d d x [\{d,<pf + \{d^f + ±(r P 2 , + r 2 P 2>2 ) + 

+ ^(/o,oP4,0 + / 2 ,oP4,2 + fifiPi,i) , (2.2) 



where 



4 

/o,o = 777(15^0 + 8w + 36w ), 
00 
4 

/2,o = g(5w - 4w - 3w ), 

/ 4)0 = 6u> - u - v Q . (2.3) 



If r 2 = f2,o = / 4 ,o — we obtain the 0(5)-invariant Hamiltonian. 

The RG dimensions of the relevant quadratic perturbations P 2i o and P 2i2 are 
respectively 21 ' 4,22,23 y 2i0 = 1/V = 1.31(1) (u is the 0(5) critical exponent associated with the 
correlation length) and y 2)2 = 1.83(1). The corresponding crossover exponents are 2jO = 1 
and <p2,2 — V2,2 V — 1-40(2). The corresponding Hamiltonian parameters and must be 
tuned to approach the MCP. The spin-0 quartic perturbation P 4 determines the leading 
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scaling corrections in the 0(5) vector model; its RG dimension is y 4) o = —0.79(2). 4 The 
spin-2 quartic perturbation P 42 turns out to be irrelevant: y^ 2 = —0.441(13). 4 The 0(5) 
FP is also unstable with respect to the perturbation P 4j4 (more generally, the O(N) FP is 
unstable against the spin-4 quartic perturbation for N > N c with N c < 3, see Refs. 24-27). 
The corresponding RG dimension y 4i4 is given by y 4i4 = ^e + 0(e 2 ) close to four dimensions. 
A naive extrapolation of this result to three dimensions, achieved by setting e = 1, suggests 
a very small value. On the other hand, 27 ' 4,23 a six-loop analysis of the 3-d perturbative series 
in the massive zero-momentum scheme gives y 4j4 = 0.189(10), while the five-loop e expan- 
sion gives y 4i4 = 0.198(11). These estimates are still moderately small, though much larger 
than the estimate y 4i4 = 1/13 obtained by naively setting e = 1 in the one-loop e-expansion 
result. 



B. Monte Carlo determination of the critical exponent t/ 4 4 

In the following we present a MC calculation of the RG dimension y 4j4 and of the crossover 
exponent 4j4 = z^/ 4)4 . As we shall see, the result will provide conclusive evidence for the 
fact that the 3-d exponent 4j4 is actually much larger that the O(e) result, in agreement 
with the high-order FT analyses. 

We consider the 5-vector model on a simple cubic lattice of size L 3 with Hamiltonian 

H = -P'£s x -s y , (2.4) 

where s x is a five-component unit vector, and the summation is extended over all nearest- 
neighbor pairs (xy). The RG dimension y 4;4 can be obtained by studying the scaling behavior 
of the cubic-symmetric perturbation 

^c = EIXi (2.5) 

x j=l 

at the critical point. Indeed, the operator P c is a particular combination of the spin-4 and 
spin-0 operators 4 and thus its scaling behavior allows the determination of y 4j4 (assuming of 
course that y 4i0 < 1/4,4)- We follow closely Ref. 26 where y 4)4 was computed in the N- vector 
model for N = 2, 3, and 4. We define 

V 5 /If 4 

M = Y,s x , R=^f, (2-6) 
and compute the correlation 

D R =(RP C )-(R)(P C ) . (2.7) 

The RG dimension |/ 4 4 is obtained from the finite-size scaling behavior of Dr at the critical 
point, since 

D R (L,P = p c )~L^ (2.8) 
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for L — > oo. 

In our simulations we used combinations of wall-cluster, 28 single-cluster, 29 and overrelax- 
ation updates. In order to determine the critical point and the standard critical exponents, 
we considered lattices of size L < 128 for several values of (3 close to 1.1813. We measured 
the Binder cumulant 

U 3 = (2.9) 

(M 2 ) 2 

and the ratio Z a /Z p , where Z p is the partition function with periodic boundary conditions 
and Z a is the partition function with periodic boundary conditions in two directions and 
antiperiodic boundary conditions in the third one (for a detailed discussion of this quantity 
see Refs. 30, 18, 17 and references therein). 

The critical value f3 c is estimated by employing the standard crossing method for the 
Binder cumulant and for Z a /Z p . We obtain 

#. = 1.18138(3), (2.10) 

where the quoted uncertainty includes both statistical and systematical errors. It is compat- 
ible with the estimate f3 c = 1.18127(10) obtained using only lattices with L < 32, indicating 
that scaling corrections are quite small. The result (2.10) is consistent with the estimate 31 
obtained from the analysis of the 21st-order high-temperature expansion of the magnetic 
susceptibility computed in Ref. 32, but much more precise. The analysis also provides the 
critical values of U 4 and of Z a /Z p : V\ = 1.069(1) and (Z a /Z p )* = 0.071(1). 

We estimate the exponent v from the L-dependence of the slope of Z a /Z p and of U4 
along the line in the (j3,L) plane on which Z a /Z p = 0.071 (which means that, for each L, 
we take the two quantities at the value of j3 where Z a /Z p takes the value 0.071). Fitting 
the slope of Z a /Z p for data with 16 < L < 128 with a simple power-law Ansatz, we obtain 
v = 0.7782(14) with x 2 /d.o.f.= 0.49, where d.o.f. indicates the number of degrees of freedom 
of the fit. If we use the data with 8 < L < 64 we obtain instead v = 0.7766(5) with 
X 2 /d.o.f.= 0.95. The systematic error due to scaling corrections should not be (much) larger 
than the difference between these two results. Roughly, assuming that scaling corrections 
are dominated by the leading one with exponent cu = — 2/4,0 = 0.79(2), the systematic error 
on v = 0.7782(15) should be given approximately by (2 U — 1) _1 times the difference, i.e. 
it should be approximately 0.002. Fitting the slope of the Binder cumulant, we obtain 
v = 0.7809(32) (resp. v = 0.7802(13)) from data with 16 < L < 128 (resp. 8 < L < 64). 
These results are consistent with those from the slope of Z a /Z p but less precise. As our final 
estimate we give 

v = 0.779(3), (2.11) 

where the error should also include the systematic uncertainty due to scaling corrections. 
This estimate clearly rules out the MC result v = 0.728(18) of Ref. 9. It is also slightly 
larger than the FT estimate 22 v = 0.762(7). Moreover, by fitting the magnetic susceptibility 
X at Z a /Z p = 0.071 we find 

77 = 0.034(1), (2.12) 



7 



TABLE I. Estimates of Dr at (3 = 1.18138. The number reported in parentheses is the 
statistical error, while the number in brackets gives the error due to the uncertainty on the estimate 

of A, 



L 


Dr 


4 


0.023979(6) [4] 


5 


0.025668(8) [6] 


6 


0.026950(10) [9] 


7 


0.028019(13) [10] 


8 


0.028942(16) 12 

V / L J 


9 


0.029733(18) [15] 


10 


0.030501(22) [18] 


11 


0.031100(25) [22] 


12 


0.031752(28) [25] 


14 


0.032900(36) [31] 


16 


0.033825(44) [38] 


18 


0.034769(53) [47] 


20 


0.035638(81) [47] 


22 


0.036461(72) [60] 


24 


0.037134(88) [70] 



where again the quoted error includes possible systematic errors due to scaling corrections. 
Performing this analysis we followed closely Refs. 30, 18, 17, where the models with 2, 3, 
and 4 components were studied. 

The determination of y±^ through Eq. (2.8) requires a very large statistics, which limits 
us to lattices with L < 24. Typically, approximately 10 9 measurements were performed 
for each lattice size. Most of the simulations were performed at f3 — 1.18127. In order 
to extrapolate the data to other values of (3 and, in particular, to /3 C , we also performed 
simulations at f3 — 1.18 and (3 = 1.183, obtaining an estimate of the derivative of Dr with 
respect to (3. In Table I we summarize the numerical results for Dr at (3 — 1.18138 ~ (3 C . 

First, we fit these data with the simple power-law ansatz (2.8). The results are sum- 
marized in Table II. The fits are stable starting from L min = 8. Also the % 2 /d.o.f. for fits 
with L min > 8 is smaller than one. In order to check the dependence of the result for 
on the numerical estimate of /3 C , we repeat these fits using the values of Dr at f3 — 1.18135 
and /3 = 1.18141, which correspond to adding or subtracting one error bar to f3 c . For ex- 
ample, for L min = 11 we obtain the estimates y 4i4 = 0.2256(19) and y 4i4 = 0.2283(19) for 
(3 = 1.18135 and (3 = 1.18141, respectively. Optimistically, we could quote the final result 
1/4,4 = 0.227(3), where the error bar includes the statistical uncertainty as well as the error 
due to the uncertainty in the estimate of f3 c . This estimate would be reliable if scaling correc- 
tions, and in particular the leading ones associated with the exponent uj = — y^o = 0.79(2), 
are suppressed. This seems to be supported by the above-reported analysis to determine (3 C 
and the critical exponents, where no evidence for the leading scaling corrections was found 
within our statistical errors. 
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TABLE II. Fits of Dr((3 = 1.18138 « (3 C ) on lattices of size L > L Ui \ n with the power-law 
Ansatz Dr = cL y4 " 4 . In the final column the \ 2 P er degree of freedom (d.o.f.) is reported. 



min C X /d.O.f. 

6 0.2333(3) 0.01778(2) 8.30 

7 0.2290(8) 0.01796(3) 2.20 

8 0.2268(10) 0.01807(4) 0.93 

9 0.2261(12) 0.01810(6) 0.92 

10 0.2252(15) 0.01815(7) 0.90 

11 0.2270(19) 0.01806(9) 0.60 



In order to be conservative, we also analyze the data allowing for a nonvanishing scaling 
correction, i.e. we fitted our data with the ansatz 

D R (L,(3 = f3 c ) = cV>^ (1 + c'L-) , (2.13) 

where we set iv = 0.79. 4 Note that the results for y 4;4 depend very little on the precise value 
for iv. The results of the fits are summarized in Table III. Only for L m[n < 8 do we get 
an amplitude d that is different from zero within error bars. However, for these fits also 
X 2 /d.o.f. is large. This suggests that for the small lattice sizes L < 8 actually subleading 
scaling corrections dominate. The results for the exponent y 4j4 obtained from L min > 7 
are completely consistent with the results reported above for the power-law ansatz without 
corrections. However, here the error bars are considerably larger. Since we are not able to 
exclude leading corrections to scaling on a firm basis, we take the error resulting from the 
fits that include leading corrections to scaling as our final estimate. Taking into account 
also the uncertainty in our estimate of (5 C we arrive at the final estimate 

y 4A = 0.23(2). (2.14) 

Note that the reported error estimate is probably quite conservative. Using the estimate 
(2.11) of v, we obtain 4 4 = y^u = 0.180(15) for the corresponding crossover exponent. 

III. RENORMALIZATION-GROUP FLOW 

The presence of a stable FP does not imply that the MCP should always correspond 
to a second-order transition controlled by the stable FP. Indeed, this happens only if the 
system is in the attraction domain of the stable FP. It is thus of interest to determine the 
attraction domain of the DFP. For this purpose we study the RG flow of the system, i.e., 
the RG trajectories along which the quartic Hamiltonian parameters u , v , and w are kept 
fixed. 

We consider the flow directly in three dimensions, using the minimal-subtraction (MS) 
scheme without e expansion. 20 In this scheme one considers the massless (critical) theory in 
dimensional regularization within the MS scheme. 33 RG functions are obtained in terms of 
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TABLE III. Fits of D R (f3 
The exponent co was fixed: oj 
reported. 



= 1.18138 « (3 C ) on lattices of size L > L m \ n with the Ansatz (2.13). 
= 0.79. In the final column the \ 2 P er degree of freedom (d.o.f.) is 



L m in c d_ X 2 /d-o-f. 

5 0.189(3) 0.0208(2) -0.321(16) 3.02 

6 0.199(4) 0.0201(3) -0.244(27) 1.87 

7 0.210(5) 0.0193(3) -0.152(37) 1.24 

8 0.221(8) 0.0185(6) -0.054(70) 0.98 

9 0.226(10) 0.0181(7) -0.001(95) 1.06 

10 0.245(15) 0.0168(10) 0.20(16) 0.83 

11 0.226(18) 0.0182(13) -0.01(20) 0.71 



the renormalized couplings u, v, w and of e = 4 — d. Subsequently e is set to its physical 
value e = 1, providing a 3-d scheme in which the 3-d RG functions are expanded in powers 
of the MS renormalized quartic couplings u, v,w. This scheme differs from the standard e 
expansion 34 in which one expands the RG functions in powers of e. The MS (3 functions 
have been computed to five loops for generic 0(ni)©0(ri2) LGW theories. 4 In the appendix 
we report the five-loop series that are used in this section to determine the RG flow of the 
0(3)00(2) case. 

The RG trajectories are obtained by solving the differential equations 

-A^ = /3>(A)XA)XA)], 
dv 



-A— = /3>(A), V (A)XA)], 
-A^=/3>(A)XA)XA)], 



(3.1) 



with A G [0, oo) and the initial conditions 



«(0) 
du 
d\ 



v(0) 



= u 0, ~T\ 



A=0 



w(0) = 0, 
dv 
dX 



A=0 



dw 
dX 



Wq. 



(3.2) 



A=0 



Note that the trajectories do not depend on the Hamiltonian parameters individually, but 
only through their dimensionless ratios. For instance, by rescaling A — > X/w , the initial 
conditions depend only on u /w and v /w . The necessary resummation of the series is 
perfomed by employing the Pade-Borel method; see, e.g., Refs. 35, 24. 

In this RG scheme the coordinates of the 0(5) FP are u* = v * = 3w* ~ 0.62. Since the 
0(5) FP has only one instability direction in the u, v, w space, RG trajectories flowing to 
the 0(5) FP separate the flow into two different domains. In order to characterize these 
domains, we consider the dimensionless ratios f2,o/ fo,o, /4,o//o,0) cf. Eq. (2.3), and A //q 
where A = uqVq — 9w^. In Figs. 3 and 4 we report the initial conditions that correspond to 
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-1.0 -0.5 0.0 0.5 1.0 

FIG. 3. Initial conditions for trajectories flowing to the 0(5) FP in the /2,o//o,o> /4,o//o,o plane. 

trajectories flowing to the 0(5) FP, respectively in terms of / 2 ,o//o,o an d /4,o//o,o an d i n terms 
of A //o and /4,o//o,o- The results reported in the figures were obtained by resumming the 
(3 functions with a [5/1] Pade-Borel approximant. We used several different values of the 
Borel-Leroy parameter b: the results are essentially independent of this parameter. 

Since the spin-4 perturbation is relevant, asymptotically the RG trajectories reaching 
the 0(5) FP lie in the plane f± = 0. One may then naively expect that the initial conditions 
for trajectories flowing to the 0(5) FP have /^o/Zo.o small. As it can be seen from Fig. 3 
this is not really true, although |/4,o//o,o| ^ |/2,o//o,o|- We also report the results in terms 
of A //oo f° r two reasons. First, A = u vo — 9u>q determines the nature (bicritical or 
tetracritical) of the phase diagram. Second, in the mean-field approximation Ao = is 
required for obtaining an 0(5)-invariant transition. Indeed, if one ignores fluctuations and 
thus neglects the kinetic term and rescales the fields as <\> — > {u/uq} 1 ^^ and ip — > (u/voY^ip, 
one obtains an 0(5)-invariant quartic potential only for Ao = 0. Note that, as show by Fig. 4, 
this mean-field condition is not actually satisfied when fluctuations are taking into account, 
although A //o is small whenever /4,o//o,o is a lso small. 

The curve reported in Fig. 3, of equation ^o/Zo.o = <7(/2,o//o,o)> separates the space of 
initial conditions into two different domains. We determine now which of them contains 
the attraction domain of the DFP, whose coordinates are u* ~ 0.688, v* « 0.745, w* = 0, 
i.e., /J ~ 1.86, /| ~ 0.204, fl « —1.43 (here, fa are the renormalized couplings related 
the Hamiltonian couplings /j,o)- A numerical analysis indicates that the DFP can only be 
reached for initial conditions that are inside the curve reported in Fig. 3, i.e., for values 
satisfying /4,o//o,o < g(f 2,0/ fo,o)- Points that satisfy the opposite inequality correspond to 
trajectories running away to infinity. These results can be rephrased in terms of Ao. A 
simple calculation shows that, if /4,o//o,o < fl , (/2,o//o,o) ) then A > 0. In other words, 
bicritical systems correspond to trajectories running away to infinity: in this case the MCP 
corresponds to a first-order transition. On the other hand, tetracritical systems may show, 
depending on the values of the parameters, a second-order multicritical behavior controlled 
by the DFP or a first-order transition. Note that the first-order nature of the bicritical point 
should have been guessed a priori: indeed, by its very nature, the DFP can only correspond 
to a tetracritical phase diagram and thus cannot be relevant for bicritical systems. 
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FIG. 4. Initial conditions for trajectories flowing to the 0(5) FP in the /4,o//o,o> &o/fo o plane. 
Points belonging to the lower dashed curve have /2,o//o,o > 0, those belonging to the upper 
continuous one have /2,o//o,o < 0. 

IV. CONCLUSIONS 

In conclusion, our MC simulations confirm that the 0(5) FP is unstable in the 3-d 
0(3) ©0(2) LGW theory. The instability of the 0(5) FP is related to the presence of the 
spin-4 quartic perturbation in the 0(3) ©0(2) LGW Hamiltonian. We obtain an accurate 
estimate of its RG dimension, 3/4,4 = 0.23(2), and of the corresponding crossover exponent, 
4 , 4 = y\^v = 0.180(15). These results are slightly larger than, but in substantial agreement 
with, the FT results. 4 The only stable FP of the 0(3) ©0(2) LGW theory is the DFP. 16,4 ' 19 
Therefore, if the system is in its attraction domain, a decoupled critical behavior should be 
observed at the MCP, with a tetracritical phase diagram. 

As argued in Refs. 8, 14, sufficiently close to the MCP, the 0(3) ©0(2) LGW Hamiltonian 
is the effective theory of a realistic SO (5) theory of high-T c superconductors, such as the 
projected SO(5) model. 6,11 The calculations of Ref. 11 indicate that realistic models have 
a bicritical phase diagram with A < 0. Therefore, the analysis that we have reported 
indicates that the MCP should correspond to a first-order transition and the phase diagram 
should be like that reported in Fig. 5. Note that first-order transitions should also be 
observed along the 0(2) and 0(3) lines close to the MCP. Therefore, the 0(2) and 0(3) 
lines in the temperature-doping phase diagram are expected to present, starting from the 
MCP, first-order transitions up to a tricritical point with an essentially mean-field critical 
behavior (apart from logarithms), and then second-order transitions in the 0(2) and 0(3) 
universality class, respectively. 

As put forward in Refs. 16, 15, 14, an effective 0(5) multicritical behavior may still be 
observed in a preasymptotic region when the effective breaking of the 0(5) symmetry, and 
in particular its spin-4 component, is small. The observation was essentially based on the 
value of the crossover exponent suggested by an 0(e) calculation 3 , i.e. 04,4 = e/26, which 
gives 04,4 ~ 0.038 when one sets e = 1. Using this estimate of 04,4, the authors of Refs. 15, 
14 argued that significant deviations from the symmetric 0(5) multicritical behavior cannot 
be observed in experiments and realistic MC simulations, since the crossover to the eventual 
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FIG. 5. Phase diagram with a first-order bicritical point. Here, T is the temperature and g a 
second relevant parameter. The thick lines correspond to first-order transitions. 

asymptotic behavior is expected only when the reduced temperature is t C ross ~ 10~ 10 , 10~ n . 
On the other hand, as shown in this paper, the actual 3-d value of 04,4 is much larger. 
Taking for granted the arguments leading to the estimate t croS s ~ 1CT 10 , one can easily see 
that, using the actual value of 04,4, the crossover reduced temperature t croS s would change 
from 10~ 10 to 36 tcross ~ 10~ 2 . Taking into account that this is a very rough estimate, which 
may easily miss one order of magnitude, we conclude that the actual 3-d value 04,4 ~ 0.18 
is sufficiently large to give rise to observable crossover effects towards the real asymptotic 
behavior even in systems with a moderately small breaking of the 0(5) symmetry, such as 
the projected S0(5) model. 11 These crossover effects should be observable without the need 
of reaching extremely small values of the reduced temperature, as argued in Refs. 14, 15. 

We would also like to stress that the value 04,4 ~ 0.18 is of the same order of magnitude 
of the crossover exponent appearing in other physical systems. For instance, in randomly 
dilute uniaxial magnetic materials, see, e.g., Refs. 24, 37, the pure Ising FP is unstable with a 
crossover exponent ps 0.11, which is substantially smaller than 4 , 4 at the 0(5) FP. In these 
systems the asymptotic critical behavior has been precisely observed both numerically 38 ' 39 
and experimentally, 40 and sizeable crossover effects from the Ising to the random-exchange 
critical behavior have been observed also in the case of small dilution, without the need of 
reaching extremely small reduced-temperature values. 

Since we have computed the initial conditions that guarantee the flow to reach the 0(5) 
FP, we can try to understand more quantitatively the behavior of the model discussed in 
Ref. 11. Using their notations, at the classical level the projected SO(5) model is equivalent 
to the LGW theory with 




5(r? 2 - 1) 
2(2 + 3t? 2 )' 




(4.1) 
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where r] is related to the different mobilities of the holes and of the magnons. The model 
is SO (5) invariant for 77 = 1 as expected. As rj decreases (a realistic value would be 15 rf = 
0.225), /2,o//o,o becomes negative and increases in absolute value. Comparing with Fig. 3, 
we see that the SO (5) breaking effects increase. An estimate of the crossover temperature 
Across can be obtained by requiring 41 5t~^ 4A ~ 1, where S somehow measures the distance 
between the initial condition and the line corresponding to the 0(5) invariant models. In 
this case, we can take 5 = |/2,o//o,o| an d thus estimate t cross ~ 1CT 1 for the realistic case 15 
rf = 0.225. Note that t cross decreases rapidly with increasing r\ and indeed one obtains 
Across ~ 10~ 3 for the model studied numerically (rf = 1/2) in Ref. 15, thereby explaining the 
apparent 0(5) behavior that was observed. Inclusion of the quantum corrections has the 
effect of increasing the SO(5) breaking. For instance, for rf = 0.225 we have 

^ = -0.724 - 4.36k, ^ = 1.26k, (4.2) 

Jo,o Jo,o 

where k = VDJI S (see Ref. 11 for the definitions) is a positive quantity. Thus, for the 
realistic case, the initial conditions are pushed even further from the line corresponding to 
the 0(5) invariant models. 



APPENDIX A: THE FIVE-LOOP SERIES OF THE (3 FUNCTIONS IN THE MS 

SCHEME 



The expansion of the (3 functions of generic 0(ni)ffiO(n 2 ) LGW theories in the modified 
minimal-subtraction (MS) renormalization scheme is known to five loops. 4 In Ref. 4 only 
the e expansions of some exponents were reported. Here we report the five-loop series of the 
(3 functions of the 0(3)©0(2) LGW theory (1.1) in powers of the MS renormalized quartic 
couplings. 

In order to renormalize the 0(3)©0(2) LGW theory (1.1) in the MS scheme, one 
sets = Z l J 2 ct) r) ip = Z 9 (u,v,w) 1/2 tp r , u = A d fi e Z u (u,v,w), v = A d ^ € Z v (u,v,w), 
wo = A d Li e Z w (u,v,w), where u,v,w are the MS renormalized quartic couplings. The renor- 
malization functions and Z UjVjW are normalized so that Z^^ m 1 and Z UjVjW rs u,v,w 
at tree level. Here A d is a <i-dependent constant given by A d = 2 d - 1 n d / 2 T(d/2). The MS (3 
functions /3 UjVjW are obtained by differentiating the renormalized couplings with respect to 
the scale //, keeping the bare couplings u , v , w fixed. The 0(3)©0(2) (3 functions are 
given by 



(3 U 


= (d 


(3 V 


= (d 


(3 W 


= (d 



l)u + ^ru 2 + 3w 2 - ^u 3 - \uw 2 -6w 3 + J2 b^ k uVw k , 

ijk 

l)v + -v 2 + -w 2 - -v 3 - —vw 2 - 9w 3 + V b'flu i v j w k , 
' 3 2 3 4 u 

l)w H — uw H — vw + 2w 2 — —u 2 w uw 2 

' 6 3 72 2 

_l w 2 w _ W _ 21^3 + y> b^uVw k . 
18 8 ^rrf %3k 

ijk 



(Al) 
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The coefficients b\j£ up to five loops, i.e. for A<i + j + k<6, are reported in Tables IV. 
In order to save space, we report them numerically, although we have their exact expressions 
in terms of fractions and of ( functions with integer argument. 
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TABLE IV. Coefficients of the (3 functions, cf. Eq. (Al). 



i,3,k 



b 



TOT 

3,1,0 
3,0,1 
2,2,0 
2,1,1 
2,0,2 
1,3,0 
1,2,1 
1,1,2 
1,0,3 
0,4,0 
0,3,1 
0,2,2 
0,1,3 
0,04 

4,1,0 

4,0,1 

3,2,0 

3,1,1 

3,0,2 

2,3,0 

2,2,1 

2,1,2 

2,0,3 

1,4,0 

1,3,1 

1,2,2 

1,1,3 

1,0,4 

0,5,0 

0,4,1 

0,3,2 

0,2,3 

0,1,4 

0,0,5 

"B^IT 

5,1,0 

5,0,1 

4,2,0 

4,1,1 

4,0,2 

3,3,0 

3,2,1 

3,1,2 

3,0,3 

2,4,0 

2,3,1 

2,2,2 

2,1,3 

2,0,4 

1,5,0 

1,4,1 

1,3,2 

1,2,3 

1,1,4 

1,0,5 

0,6,0 

0,5,1 

0,4,2 

0,3,3 

0,2,4 

0,1,5 

0,0,6 



(u) 
ijk 



b 



h.95643 





2.21875 



-1.08333 
38.6747 



0.125 

8 

13.9123 
-27.3529 





-5.56225 



1.18623 
-188.898 



-0.144409 
-19.1406 
-173.188 





-0.489598 
-10.9914 
-91.5308 
-120.043 
156.207 









6.47354 



2.43996 
1146.82 





1.98376 

35.922 

1294.73 





-0.142137 
34.0863 
659.256 
1776.58 





1.27762 
37.2454 

257.6 
943.812 
1108.93 



0) 
ijk_ 



b 



0" 






0.234375 



-2.03125 

15 
4.99347 


3.46875 
54.262 
20.4465 











-1.0098 



-0.0839266 
-21.3002 



1.81816 
-34.1777 
-169.835 
-21.9072 


-7.99393 
-247.731 
- 239.138 
-191.964 











2.81164 



-0.068611 
77.9976 



3.09444 
57.0315 
518.911 





4.13976 
66.1707 
1136.61 
1751.3 
120.141 


10.0207 
1426.49 
1657.51 
2594.22 
1879.79 



(w) 
ijk 



o 

0.902778 



4.44611 



I. 38889 
14.0051 



0.652778 
3.50134 

II. 6624 
12.8874 




-2.91484 



-17.0442 



-2.56238 
-51.1747 





-2.57945 
-41.1293 
-108.727 


-1.99192 
-12.3167 
-39.1403 
-94.7246 
-100.117 





12.3408 



80.795 



11.0339 
249.36 





5.96446 
124.802 
617.85 





9.84045 
119.023 
636.697 
1155.64 


7.99517 
54.8793 
178.161 
484.147 
1031.67 
852.872 
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